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We investigate the problem of dynamical gap generation in suspended graphene by long-range 
Coulomb interactions at strong coupling with Dyson-Schwinger equations. Including renormal- 
O ization effects on the Fermi velocity we obtain a critical coupling constant a c = 2.85 which is 

larger than the bare coupling Oo = 2.19 of suspended graphene. This suggests that at low ener- 
gies the running of the Fermi velocity prevents the emergence of a gapped phase. Our calculation 
is thus in agreement with the experimental observation that suspended graphene remains in the 
g semimetal phase. 
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1. Introduction 

Due to its unusual electronic properties, graphene, a monolayer of carbon atoms on a two- 
dimensional honeycomb lattice [1], has been intensely studied both theoretically and experimen- 
tally over the last years (see, for example, Refs. [2, 3, 4] for overviews). In particular, much work 
has been devoted to understanding the effects of electron-electron interactions, both in the contin- 
uum, mainly via Dyson-Schwinger and renormalization group equations [5, 6, 7, 8, 9, 10, 11, 12, 
13, 14, 15], and on the lattice [16, 17, 18, 19, 20]. The majority of the theoretical efforts were 
thereby focussed on the question whether the long-range Coulomb interactions might trigger the 
dynamical generation of a gap in graphene at strong coupling similar to chiral symmetry break- 
ing in QCD. This possibility is motivated by the potentially large value of the effective coupling 
a = e 2 /ehvF, where vf ~ 10 6 m/s c/300 is the Fermi velocity of graphene on a substrate with 
dielectric constant e. If the Fermi velocity of suspended graphene with e = 1 is of the same order, 
this could lead to the formation of a condensate corresponding to a transition into an insulating 
phase above some critical coupling a c . 

In a series of theoretical works, various authors have calculated the value of the critical cou- 
pling that would correspond to the semimetal-insulator transition: Gamayun et al. found a c = 0.92 
by solving a simplified gap equation in which radiative corrections to the Fermi velocity are ne- 
glected [10]; in [12], an approximation for the photon polarization has been used, leading to the 
critical value a c = 1.13; renormalization group calculations at two loops yield a c = 0.833 [21]; on 
the lattice, calculations performed using a 'standard' square lattice found the value a c = 1 .08 ± 0.05 
[17], whereas simulations carried out directly on a hexagonal (physical) lattice yield a c = 0.9 ± 0.2 
[19]. All these values lie below the bare coupling constant of suspended graphene, Ofo = 2.19. 

On the other hand, the experimentally measured conductivity shows that suspended graphene 
is not in the insulating phase [4, 22]. It is observed, however, that the strong Coulomb interactions 
induce a charge-carrier density dependent renormalization of the Fermi velocity corresponding to 
a reshaping of the Dirac cones. This renormalization peaks at half-filling where it leads to an 
increase by about a factor of three in the Fermi velocity and thus a corresponding decrease in 
the renormalized effective coupling. A strong renormalization of Vf was predicted already in an 
early perturbative RG study [5] where it was concluded that the Coulomb interactions would in 
fact induce a logarithmic increase in the Fermi velocity without bound until retardation effects 
eventually become important and the instantaneous Coulomb approximation breaks down. 

While lattice simulations can yield reliable results in strongly coupled theories, extrapolations 
to the infinite volume and chiral limits are often a numerical challenge. This is also very much true 
for the simulations of the electronic properties of graphene. Therefore, functional methods in which 
these limits are accessed more easily can provide important alternative nonperturbative formula- 
tions. Moreover, the inevitable approximations in these formulations can be validated from simula- 
tions with sufficiently large masses in finite volumes. Here we employ the Dyson-Schwinger equa- 
tions. Within this approach, substantial progress has been made over the past decades, both in in 
QCD and QED3 (see, for example, Refs. [23, 24, 25, 26, 27] and references therein). Extending 
the study of Ref. [10], we first calculate the running Fermi velocity in the bare vertex approxima- 
tion and with a one-loop (frequency dependent) photon polarization. In a second step, we solve 
the gap equation at the critical point by explicitly including the dressing function corresponding to 
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this running Fermi velocity. For the critical coupling of the semimetal-insulator transition we then 
obtain a c = 2.85 [28] which is larger than the bare coupling of suspended graphene 0!o = 2.19, 
indicating the persistence of the semimetal phase. Our results thus agree with the experimental 
findings [22]. Finally, in order to compare our results with the lattice simulations which must in- 
clude a sufficiently large mass term in a finite volume, we have also investigated the dependence of 
the corresponding pseudocritical coupling on such a fermion mass term which acts as an explicitly 
symmetry breaking external field as the current quark masses do in QCD. 

2. Running Fermi velocity and critical point analysis 

We start by briefly sketching the details of the continuum model of graphene employed in this 
work. At the level of the Lagrangian, including the Coulomb interaction corresponds to accommo- 
dating the boson field that propagates in (3+1) dimensions into the (2+1) dimensional Lagrangian of 
the massless Dirac fermions. This is done via dimensional reduction, which gives rise to a so-called 
reduced gauge theory (we refer to Refs. [9, 10, 28] for a complete description and derivation). The 
dimensionally reduced action of 2D Dirac fermions interacting through a three-dimensional gauge 
field is given by 

y = J dtd 2 r\j/(t,r) [ifd t - ivffdi] y/(t,r) 

~hl dtdt ' d2Td27 ' 5 ( ? ~ ?/ ) y (* . ?) y° y (t, r) ■ _, 1 - . xjr (t' , r ') / Y (?' ,f' ) . (2.1) 

Here, the spatial index is i = 1,2, and s = 1, -;Nf is a flavor quantum number with Nf = 2 for 
monolayer graphene referring to the physical spin degree of freedom. denotes the fermion 

(conjugate) field, vj? is the Fermi velocity, and e is the dielectric constant of the medium (for 
suspended graphene 6 = 1, whereas on top of a substrate one uses e = (1 + k)/2 typically with 
K ps 4 for Si0 2 or K « 10 for SiC). 

From the action in Eq. (2.1) one obtains the fermion Dyson-Schwinger (gap) equation via 
standard functional derivation. In Minkowski space we have 

S- l (po,p) =S l (p ,p)-ief J ^S(ko,k)r(k,-p)D(po-k ,p-k), (2.2) 

where F(k, —p) is the fully dressed fermion-photon vertex and the fermion propagator is given by 
(without wavefunction renormalization, details will be provided elsewhere [28]): 

«/ ,_ fpo-VFA{p)fpj-k{p) 
" [P) pl-v 2 F A*(p)pi-A*(p) ' 

where A is the Fermi velocity dressing function, and A is the fermion mass function. The Coulomb 
propagator in the random-phase approximation (RPA), 

%,j)= M l! - v (2-4) 
\q\ +U(q Q ,q) 

includes polarization effects with the one-loop polarization function [5], 

U(q Q ,q) = ^—J^. (2.5) 
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Figure 1: The Fermi velocity dressing function A in GWA with RPA Coulomb interactions, A ~ %/a, 
Nf = 2 for monolayer graphene, and a — 0.8, 2.0 and 2.85 to represent the presence of a substrate, suspended 
graphene, and criticality for the semimetal-insulator transition, respectively, from bottom to top. 



Following the usual procedure, we recast the Dyson-Schwinger equation Eq. (2.2) into a set of 
coupled integral equations for the dressing functions A(p) and A(p). Neglecting a possible weak 
energy dependence in A and A, and in a bare vertex truncation F — > F^ = ey°, one obtains [10]: 



A(p) = 1 + 



do) 
2% 



d 2 k p-k 



A(k) 



(2n) 2 p 2 (o 2 + v 2 A 2 (k)k 2 +A 2 (k) 



—D((D,p-k), 



e 2 r da 



d 2 k 



A(k) 



(2k) 2 co 2 + v 2 A 2 (k)k 2 + A 2 (k) 



D(co,p-k). 



(2.6) 



(2.7) 



In order to evaluate these equations at the critical coupling we apply bifurcation theory [29], which 
amounts to linearizing Eqs. (2.6, 2.7) around the point at which the nontrivial solution for the mass 
function bifurcates away from the trivial one. In this case the integral equation for A decouples. 
Furthermore employing a GW approximation as in Ref. [30] amounts to setting A = 1 on the 
r.h.s. of Eq. (2.6) such that the corresponding equation is reduced to solving a single integral. The 
integration over the (Euclidean) frequency co can then be performed analytically [10], yielding 
(now with p = \p\): 



A(p) = l + 



v F e p 2 



d k p ■ k 
{2x) 2 k\p-k 



(2.8) 



where g = e 2 Nf7i/(4evf) and the piecewise function J(z,g) with z 



{z 2 -\)[n-gc{z)]+zg 2 c{g) 



p — k\ is given by 
' 2 arccosh (x) / Vx 2 — 1 , x > 1 



z 2 +Z 2 - 



with c(x) 



2arccos(^)/-\/l — x 2 
2 



x < 1 
x= 1 



(2.9) 



At this point, it is common to further approximate the rather complicated integration over the angle 
between p and k [10]. Here we proceed differently using z as an integration variable to encode the 
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dependence on this angle. Then the standard angular and radial integrals are replaced by 

z z 

A 1 , ' ' A 



" / dz= j 1 dz I dp + I " dz I dp, (2.10) 



JL 







[+3 

where j8 = k/p. With this trick, the integral over z in Eq. (2.8) nicely splits into intervals that 
correspond to the definition domains of the function J in Eq. (2.9). A straightforward calculation 
then leads to the following complete result for the Fermi velocity dressing function in the GW 
approximation with RPA polarization effects in the Coulomb interaction: 

A(p) = 1 + { - [n - 2g + c(g) (g 2 - 1 )] In | + f(g) + (n + 4g) ln2 - 4g} , (2.11) 

where the physical cutoff A is determined by the size of the Brillouin zone. The logarithmic term 
comes from the leading one-loop correction and was obtained in Ref. [30] already. Here we include 
all nonsingular terms which were not given therein. They can be expressed in terms of the function 

^= arccoshsln ( 2^f± ) 




— 71 arctan , _ > . (2.12) 

The sum in first term on the right, involving the coefficients 

8 n = f 1 dz(l-z 2 ) n ~ 3/2 arccosz and v„= / dtt(\ -t 2 )' 1 ^ 12 arccos-, (2.13) 
Jo Jq t 

can be evaluated numerically. The full GWA result (2.1 1) for the Fermi velocity renormalization 
function at various values of the bare effective coupling a = e 2 /evf is shown in Fig. 1. The 
logarithmic singularity at p 2 — > is rather weak at weak coupling and the nonsingular part of 
the radiative corrections remains small as well. This changes when a « 2 is used for suspended 
graphene. Both logarithmic and regular contributions are considerably larger and a correspondingly 
strong renormalization of the Fermi velocity results. If we simply use VfA{pi/vf) to estimate the 
renormalized Fermi velocity at nonzero chemical potential pi for a finite charge-carrier density this 
result describes the recent suspended-graphene experiments [22] so well that a more quantitative 
comparison maybe worthwhile [28]. The large effect also shows that a nontrivial dressing function 
A(p) should be included in the analysis of the semimetal-insulator transition at strong coupling, 
instead of using A = 1 as in [10]. We will turn to this analysis next. 

Having calculated the function A in GWA, we now use this result in the bifurcation analysis of 
the gap equation around the critical point. At leading order in A the gap equation then reduces to: 

We proceed by inserting the Fermi velocity dressing function Eq. (2.11) into the above equation, 
and solve the resulting equation by using standard numerical methods. Note that in this case, not 
only the logarithmic but also the regular contributions to A become important at strong coupling. 
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Figure 2: Mass dependence of the critical coupling for Nf = 2. 



For the critical coupling we then obtain a c = 2.85 which is larger than the bare coupling 
do = 2.19 with e = 1 for suspended graphene. This suggests that the Coulomb interactions might 
not be strong enough to generate a gapped phase in the electronic spectrum of suspended graphene. 
Our result is thus in contrast to the many previous studies mentioned in the introduction. The 
values predicted there all lie below the bare coupling for suspended graphene. In particular, our 
results are directly comparable to the work of Gamayun et al. in Ref. [10] in which the same gap 
equation with the same RPA improved Coulomb interaction was solved, however, with A = 1, i.e., 
with neglecting Fermi velocity renormalization effects, leading to a c = 0.92. 

As it turns out, lattice simulations at present also appear to underestimate the value of the 
critical coupling [17, 20]. Hence, in order to mimic the lattice setup, where a non-zero fermion 
mass term must be used in a finite volume, 1 we have introduced such an explicitly symmetry 
breaking term as a seed in the gap equation (2.14) here as well, and calculated the corresponding 
pseudocritical couplings. Analyzing the results displayed in Fig. 2, we notice that the diagram 
contains two distinct regions, each of them displaying a logarithmic dependence of the critical 
coupling on the bare mass (dashed lines in the figure). The scale separating the two regions, where 
the dynamically generated mass is bigger or smaller than the bare fermion mass, is of the order of 
magnitude mo/A = 10~ 2 . Note however that a proper comparison with the lattice requires taking 
into account the different UV scales, as graphene lattice simulations have been performed on both 
hexagonal [20] and square [17] lattices. Further studies will be necessary to clarify the issue. 

3. Summary and outlook 

In this work, we have provided an estimate of the critical coupling for the semimetal-insulator 
phase transition in graphene. We have employed the Dyson-Schwinger equation for the fermion 

'Finite- volume effects must also be studied carefully on the lattice. Here QED3 serves as an important example 
where it has been shown what a significant impact these can have, see for example [31] and the references therein. 
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propagator in the Dirac cone approximation with a bare vertex and Coulomb interactions including 
frequency dependent polarization effects at the RPA level. Importantly, however, we have explicitly 
taken into account the GWA running of the Fermi velocity. In this setting, we have obtained a criti- 
cal coupling of a c = 2.85, which is larger than the bare coupling Ofo = 2. 19 for suspended graphene. 
This allows us to conclude that the logarithmic growth of the Fermi velocity near charge neutrality 
weakens the Coulomb interaction preventing the emergence of a gapped phase, in agreement with 
the experimental observations [22]. 

Despite the qualitative agreement with experimental observations, a complete understanding 
of the physical picture calls for a full nonperturbative investigation, i.e., to solve the coupled system 
of integral equations for the photon and fermion propagators simultaneously. It is likely that includ- 
ing the nonperturbative gauge propagator and improving the vertex functions may lead to further 
corrections of the critical coupling. We recall that nonperturbative contributions to both the photon 
propagator and vertex functions lead to significant corrections in ordinary QED3, of the order of 
twenty percent for the photon propagator and ten percent for the vertex [27, 32]. These effects, 
along with finite volume and mass corrections, might in the end together explain the discrepancy 
between experiment and lattice simulations [17, 20]. 

Moreover, induced four-fermion contact interactions might become important at strong cou- 
pling as well. Such local four-fermion interaction terms as allowed by the symmetries should be 
considered, e.g., to model electron-phonon interactions, in addition to the long-range Coulomb 
fields in the effective low-energy continuum theory for the electronic excitations of graphene in 
the future. 2 Finally, another interesting problem is the inclusion a chemical potential for a finite 
density of charge carriers as induced experimentally, e.g., by chemical doping. While the fermion- 
sign problem makes this problem hard to address with lattice simulations, the modifications to 
the Dyson-Schwinger equations are quite straightforward. The resulting effects of a finite charge- 
carrier density on the running Fermi velocity are the subject of an upcoming publication. 
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